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Abstract. We present a wide class of models which realise a bounce in a spatially flat 
Friedmann universe in standard General Relativity. The key ingredient of the theories we 
consider is a noncanonical, minimally coupled scalar field belonging to the class of theories 
with Kinetic Gravity Braiding / Galileon-like self-couplings. In these models, the universe 
smoothly evolves from contraction to expansion, suffering neither from ghosts nor gradient 
instabilities around the turning point. The end-point of the evolution can be a standard 
radiation-domination era or an inflationary phase. We formulate necessary restrictions for 
Lagrangians needed to obtain a healthy bounce and illustrate our results with phase portraits 
for simple systems including the recently proposed Galilean Genesis scenario. 
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1 Introduction 

Was there a beginning of time? Or in more technical words: was there a spacelike boundary of 
the quasiclassical spacetime, beyond which the universe was necessarily in a strong quantum 
gravity regime? If there was a beginning, was the universe collapsing or expanding immedi- 
ately afterwards? Was the universe born infinitely (Planckian) small or infinitely large? If the 
universe experienced an early period of inflation, as all observations currently suggest, what 
happened before inflation [1, 2]. All these basic questions are of fundamental importance and 
remain interesting even if one disregards the possible observational consequences. 

Asking these questions has already led to unexpected discoveries. Indeed the first cos- 
mological model with a quasi-de-Sitter stage [3] and the cosmological perturbations within 
it [4] 1 was invented as an attempt to explain how the universe could have avoided the initial 
singularity. Twenty years later it was nonetheless proven that inflation with canonical kinetic 
terms did not solve the singularity problem [6]. 

Bouncing cosmological models provide an interesting possible alternative to standard 
Big-Bang cosmology e.g. [7-15]. 2 However, such models are plagued by significant obstacles 
and frequently exhibit pathological behavior, for criticism see e.g. [20]. Many of these patho- 
logies result from the need [21-23] for the energy- momentum tensor to violate the null energy 
condition ("NEC": T flu n^n L/ ^ 0, for all null vectors n^) to bounce a spatially flat Friedmann 
universe. 3 The question of whether a stable violation of the NEC is possible is also crucial 
for the understanding of the current accelerated expansion of the universe, see e.g. [26-29]. It 
has proven extremely difficult to construct local field-theoretic models violating the NEC in 
the context of standard general relativity. Until recently it was the common wisdom that any 

1 See [5] for the English translation. 

2 For some recent reviews on bouncing cosmologies see [16-19]. 

3 Note that it is not difficult to construct a non-singular bouncing universe with positive spatial curvature 
which would compensate the positive energy density at the bounce. For that one does not need to violate 
the NEC, but the universe will be substantially closed. A classic example of such a non-singular universe is 
given in [24], where one can also find a slow-roll regime in an m 2 (j) 2 potential in a contracting universe. For 
the most recent and simple example, see [25]. 
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violation of the NEC leads to internal inconsistencies such as: ghost degrees of freedom and 
gradient instabilities (i.e. imaginary speed of sound) [30-38]. 4 One may try to avoid these 
problems in the context of effective field theory by adding higher derivatives to the scalar- field 
action see e.g. [45-48]. This approach is useful for a better systematic understanding of the 
perturbative expansion but cannot elucidate the behavior of the cosmological background. 
Inclusion of generic higher derivative terms results in new ghost degrees of freedom [49], 5 
which, in some cases, may be moved above the naive cut off. This procedure is rather delic- 
ate [51-53] and implicitly incorporates the assumption that a healthy UV completion which 
takes care of the ghosts and gradient instabilities exists, which may not be the case [54]. If 
dangerous ghost-like instabilities are ignored, i.e. one assumes that some unknown mechan- 
ism exists which would cure them, it is easy to realize a bouncing scenario, for example, in 
[55-57]. 

Recently the situation has changed with the rediscovery of scalar-field theories with 
higher derivatives in the action, but which maintain second order equations of motion. These 
higher-derivative theories possess only the three standard degrees of freedom — the graviton 
and the scalar — allowing the theories to circumvent the conclusions of [49]. Originally these 
second-order theories were derived in [58]. 6 The simplest versions of these actions arise in 
certain modifications of gravity [64] when considered in the decoupling limit [65-67] and 
were then generalized to the so-called Galileons in [68] for the fixed Minkowski metric and 
in [69, 70] for dynamical spacetime. 7 It is the presence of higher derivatives in the action 
and the corresponding kinetic mixing between the scalar and the metric which allows for 
a stable violation of the NEC [72, 73]. The simplest class of these second-order theories 
which are minimally coupled i.e. which do not involve any direct couplings to the Riemann 
tensor, but still possess this kinetic mixing/braiding was introduced in [73] under the name 
of Kinetic Gravity Braiding. 8 This class of theories is also singled out from the most general 
second-order theories by the correspondence with the hydrodynamics of imperfect fluids [75] . 
Indeed it is this imperfection which allows these theories to avoid the pathologies pointed 
out in [32] in the case of perfect fluids. 

In the current paper we will use exactly this class of theories with Kinetic Gravity 
Braiding to study bouncing cosmologies. The possibility of a healthy bounce in a particular 
model of this class, the so-called Conformal Galileon, was mentioned in [76] where the authors 
concentrated on an always-expanding and superaccelerating stage of the evolution of the 
universe: Galilean Genesis. The details of this model were further investigated in [77] and 
[78], with the latter work focusing on the bouncing solutions of the Conformal Galileon. 

In this paper, we demonstrate that manifestly stable bouncing solutions in generic 
theories with Kinetic Gravity Braiding are rather common. For simplicity, we concentrate on 
two broad categories of models: In one, the shift-symmetric scalar-field evolves in the presence 

4 However, see [39-42]. Moreover, it is unclear whether ghost-like instabilities are present in scalar theories 
with constraints such as [43]. These theories can violate the NEC without any gradient instabilities [43] and 
can realise an oscillating nonsingular universe [44]. 

5 For a modern and detailed discussion see [50]. 

6 A subclass of these models was also considered later in a different context in [59, 60]. The result of [58] 
was independently rediscovered in [61]. The equivalence of these results was shown in [62]. The original 
Horndeski's theory was recalled for the first time in modern literature — "resurrected" — in [63]. 

7 Note that general relativity does not allow the theory to maintain the Galilean symmetry in curved 
spacetime in a manifestly self-consistent fashion. For a reduced notion of Galileon symmetry in curved 
spacetime, see [71]. 

8 See also Ref. [74] where this class of models was studied slightly later under the name of G-inflation. 
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of external hydrodynamical matter. In the other, the scalar field is not shift symmetric but 
it is the only source of energy density in the universe. For both categories we derive the 
generic conditions for high-frequency stability around the bounce point. We have found that 
it is not difficult to avoid ghosts and gradient instabilities around the bounce. However, we 
would like to stress that we have not studied the issues related to a possible strong coupling 
of the scalar perturbations. To illustrate our general analysis we study the phase portraits 
of particular systems. 

We find that in most of the considered models, even though the bounces are healthy, 
the trajectories are not free of problems. Typically, pressure singularities [79] or Big Rips 
[80] 9 are present at one of the ends of the trajectories, either in the past or in the future 
with respect to the bounce. At best, the trajectories evolve to or from regions of phase space 
where the sound speed of the perturbations becomes imaginary. This means that we cannot 
really trust the background dynamics which we calculate for the whole evolutionary history. 
As a result, our bouncing models, as they stand, do not resolve the initial cosmological 
singularity. Moreover, the amount of expansion or contraction between any such singularity 
and the bounce is rather limited. We have not touched on the generation of cosmological 
perturbations here, but it is clear that the mechanism would have to be different from the 
inflationary one. 

However, we have found a category of models (the hot G-bounce) which bounce and 
then evolve to a healthy and stable future, where the scalar has redshifted away and any 
other accompanying fluid present in the universe dominates the dynamics. As such, in these 
models we have a bounce followed by a hot Big-Bang, or possibly an inflationary period. In 
all such trajectories, at some point before the bounce, the sound speed is imaginary. On the 
other hand, the physical energy scale at which this occurs can be made much smaller than 
the Planck energy. One could hope that this problem could be resolved without recourse to 
quantum gravity, but by modifying the scalar model in some way. We would like to note 
that the presence of stable bounces is so generic that we suspect that given sufficient effort 
one should be able to construct models which remain under control over the whole history, 
providing a never-singular evolution for the early universe. 

The remainder of this paper is organized as follows: We begin in §2 by introducing 
the model and discussing the general properties of a braided scalar and its perturbations. 
In §3, we discuss bounces in the presence of external matter in theories which are shift 
symmetric. This simplifies the phase space considerably and allows us to find a prescription 
for models which bounce in a healthy fashion. To illustrate this class, we introduce the hot 
G-bounce model in §3.1, which evolves to a radiation-domination era following the bounce. 
We present a selection of other models in this class in §3.2. In §4, we present the general 
properties required to build a successful bouncing model with negligible external matter and 
then discuss in detail the dynamics of the bounces in the conformal Galileon model in §4.1. 
We conclude in §5. 

2 General Properties 

In order to realise a bounce in a spatially flat Friedmann universe, the theory has to violate 
the null energy condition ("NEC") [21]. The simplest system capable of exhibiting a large 
violation of the NEC without any linear instabilities is a kinetically braided (or galileon) 
scalar field, which we have denoted as <p. 

9 For an earlier and more detailed discussion see [81]. 
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In order to aid the reader, we recap the main equations describing background evolution 
in cosmological models with Kinetic Gravity Braiding and simple Galileons. In addition, we 
provide the formulae determining the high-frequency stability of the model, i.e. inequalities 
which need to be satisfied to prevent ghost and gradient instabilities. We present the results 
in the form of [73] and [75]. We will assume that a spatially flat Friedmann universe is 
filled with the the scalar field <f> and some hydrodynamical matter with energy density p and 
pressure p = wp. 

The gravitational part of the action is given by the standard Einstein-Hilbert term while 
the action for the scalar is 10 

S^ = J d 4 x^=5 [K {<f>, X) + G (0, X) D<j>] , (2.1) 

where 

X = ^V^V^, (2.2) 

and Vu denotes a covariant derivative, so that □ = g >MU 'V^'V u . Further it is convenient to 
introduce the diffusivity, k, which measures the deviation of the energy-momentum tensor 
from the perfect-fluid form: 

K = 2XG tX • (2.3) 

Here and throughout the paper we use the notation ( ) x = 9( )/dX. In [75], it was shown 
that <fi plays the role of an effective mass or chemical potential. We will use notation 

m = <j), (2.4) 

henceforth. The total pressure of the scalar field in the reference frame moving with the 
scalar, with velocity = d^/m, is 

V = K -m 2 G^- Km. (2.5) 

In a Friedmann universe with the Hubble parameter H, the shift charge density is given by 

n = K, m - 2mG^ + 3Hk . (2.6) 

The total energy density is given by an analogue of the thermodynamical Euler relation 

£ = ran — V — kiti = m {K )Tn — mG^) — K + 3Hmn . (2-7) 

The first Friedmann equation reads 

H 2 = l -{S + p), (2.8) 
while for the second Friedmann equation we have 

H = --(£ + p + V +p) -- (ra - ran - (p + p)) . (2.9) 
The equation of motion for the scalar field can be written in the form 

Drh + 3n (h~ ^kwA + £^ = -k( P + P ), (2.10) 



10 Unless explicitly stated otherwise, we use reduced Planck units where Mpi = (8ttGn) = 1 and the 
metric signature convention (H — ). 
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where the positivity of the quantity D implies that the perturbations of 4> are n °t ghosts, 

3 

D = n m + + -k 2 > 0. (2.11) 

In the above, we have taken the Hubble parameter H as an independent variable in the differ- 
entiation, n = n (4>,m,H) and £ = £ (4>,m,H). Finally the absence of gradient instabilities 
requires the speed of propagation of acoustic perturbations of eft to be real: 

n n + k + k(H — Km/2) n + k + k(H — Km/2) 

Dm ' S tm -3k(H- Km/2) { ' } 

The background dynamics is described by four first-order equations: m = (ft, the con- 
tinuity equation for matter p + 3H (p + p) = 0, Eq. (2.9), Eq. (2.10), plus a constraint that 
is the first Friedmann equation (2.8). The system moves on a 3d hypersurface in the phase 
space (4>,m, p, H). There are two important cases in which the dynamics greatly simplify: 

• the scalar-field Lagrangian is symmetric with respect to constant shifts in field space: 

<P — > 4> + c 

• one can neglect external matter: p = 

In both these cases, the system moves on a 2d surface. In the following analysis, we will 
concentrate on these two cases in turn. 

3 The Bounce in Shift-Symmetric Theories with External Matter 

Here, for simplicity, we are going to make an assumption that p = wp with w = const. 
The phase space for this dynamical system is a 2d surface in (m,p,H). Usually one would 
parametrize this surface by the coordinates (m,p). However, here it is not very convenient 
because the constraint (2.8) has complicated branches, and this coordinatisation works only 
locally. 11 In fact the best choice of coordinates here is (m,H). This is also very helpful for 
our purpose of analysing the bouncing solutions. Indeed it is easy to solve the Friedmann 
equation with respect to p and substitute the result into Eq. (2.9) and Eq. (2.10). We have 

p (m, H) = 3H 2 — 2>KmH — mK m + K . (3.1) 

As we assume standard stable matter such as dust, radiation, etc., the positivity of p restricts 
the region(s) in phase space to which the original physical system can evolve. Note that curves 
p (m, H) = correspond to the Id phase space describing dynamics of the shift-symmetric 
scalar field in a universe containing no accompanying matter. Thus these curves are solutions 
of the equations of motion and the trajectories never cross these boundaries. 

In shift-symmetric systems, the charge density (2.6) reduces to n (m, H) = K^ m + 3Hk 
and is conserved. Having eliminated p, we can write the dynamics of the universe as the 
first-order autonomous system: 

3k (1 + w) (K + 3H 2 - mn) - 3n(2H - Km) 

m = — — - — 2^rd — 1 ■ <3 - 2) 

3kuH + n m i{\ + w) (K + 3H 2 — mn) + nm) 

H = ^ n >- '- . (3.3) 

2iTijYi ~\~ 3k 



ii 



For the discussion of complications arising due to a similar branching, see e.g. [82] . 
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Owing to shift-charge conservation, this system is integrable i.e. possesses the following first 
integral 

I(rn,H) = — = TT9 ( K ,rn + 3HK) 1+W = 
y ' ' p SH 2 + K — m {K tTn + SHk) v ; 

Each trajectory can be parameterised by the value Iq of the first integral and is given by the 
solution of I (m, H) = Xq. Further one can substitute m (m, H) into (2.12) and obtain 

2 ( = J.M'-tf.watJ.M 

4mD z 

where 

f 2 ( m ) = 6(5 + 3w)kk m , 

fi (m) = 2k (12k 2 — 3m (1 + 2>w)KK, m + 8K mm ) , 

/o (m) = -3mK A + 6kk m ((1 + - mwK^ m ) + 4K <m K >mm + k 2 (6K >m - 2mK mm ) . 

c 2 (m, i?) = is a quadratic equation with respect to H , so that it appears always possible 
to chose k, w and K such that the sound speed is positive in the whole region where D > 0. 

At the bounce, denoted with the subscript 0, the positive energy density of the fluid is 
compensated for by the negative energy density of the scalar field: 12 

p = -£ = (K-mK :m )\ >0. (3.6) 

Further at the bounce 



\l+w 
m ) 



K - mK^ m 



(3.7) 



o 



for some Xo. 13 

In order to realise a bounce and not a recollapse one has to require that 

H = K ' mm ^ + W ) K ~ wmK )Tn ) > o ^ 8 ^ 
2K mm -\- 3k 2 

Taking into account the no ghost inequality 2Dq = 2K }jnm +3k 2 > we arrive at the following 
two options at mo : 

3 

--K 2 <K mm <0 and ((1 + w) K - wmK, m ) > , (3.9) 

if,mm>0 and ((1 + w) K - wmK >m ) < . (3.10) 



12 One could expect that this cancellation might be a substantial source of isocurvature perturbations. 
13 From here one can see that for normal matter, with w < 1, the trajectories which can bounce (or crunch) 
and go to, or appear from H — ±oo with finite m, build measure zero. Indeed, 

lim T(m,H) = 3 w k 1+w H w - 1 , 

H—>oo 

thus I (m, oo) = 3k 2 for w — l and I (m, oo) = for normal matter with w < 1. Thus for these trajectories 
To — , and (3.7) generically has a finite countable number of solutions. 
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Negative K mm at the bounce generically implies that somewhere in phase space D < 0. Note 
that the boundary of D = is a pressure singularity [73, 75, 79]. Finally we have to require 
the absence of gradient instabilities: c 2 > 0. At the bounce this implies that 

mfo(m) = — 3m 2 k 4 + GmK,K >m ((1 + w)K — mwK ym ) + (3-11) 
+ 4mlf m if mm + n 2 m (6if m - 2mK imm ) > . 

Finally we would like to stress that we are only considering the stability with respect to 
high-frequency perturbations. In particular, we have ignored the possible tachyonic masses 
and related instabilities, e.g. Jeans instability. Note that perturbations 5<f> and 5p diagonalise 
the equations for cosmological perturbations in the short-wavelength limit only. 

3.1 The Hot G-Bounce 

In this section we will present an example from a class of models which bounce in a healthy 
manner and then proceed to evolve to a phase mimicking radiation domination. We will 
model the bounce as occurring in the presence of radiation. However, as we will show, the 
presence of external matter is not actually necessary and its equation of state is largely 
insignificant. This sort of bouncing trajectory is possible even in a universe containing only 
the kinetically braided scalar and no external fluid. 

In particular, we will analyse the model which, in addition to standard Einstein grav- 
ity and an external radiation fluid with equation of state parameter w = 1/3, contains a 
minimally coupled kinetically braided scalar with the Lagrangian functions 

2 2 6 

K = X — a 2 X 3 = ™y~ and n = 2xX = xm 2 . (3.12) 

In the function K above, the term X 2 is missing, therefore it appears that the model is fine 
tuned. However, we have only picked such a model since it is the presence of the X 3 term 
that is the essential component driving the dynamics discussed below and such a simplified 
choice makes the calculations somewhat more tractable. Adding the X 2 term back in — or 
indeed even higher powers of X — does not prevent the model from having largely the same 
behaviour, provided the coefficients are appropriately chosen. The coefficient of X 2 needs 
to be somewhat smaller than a (but can be of the same order) without much change to our 
conclusions (see figure 2 for the phase portrait for a hot G-bounce model including an X 2 
term). 

The charge density (2.6) is 

n = m ^1 — -a 2 m 4 + 3Hxm^j . (3.13) 

The equations of motion can be obtained from Eqs (3.2). It is convenient to rescale time in 
these equations: 

t = r^, (3.14) 

so that the chemical potential m = <fi and the Hubble parameter H = a/a will also corres- 
pondingly rescale to become dimensionless, 



Here, we should also comment on the constant a. It has dimension M~ 4 . Since one 
would expect this mass scale, M a , to typically be below Mpi, this implies that in Planck 
units a is some number larger than one. As we will see, the bounce and associated dynamics 
occur in the region fj, ~ 1, h < 1. This means that by picking a low scale for a, we can 
make the physical scale of the bounce significantly below Mpi and thus ensure that gravity 
is under control, away from its quantum regime. 

In these new variables, the system (3.2) takes the form 

_ (15^ 4 - 4) (48/i 2 + 4/i 2 + fi 6 ) - 96/3/t/i (l2h 2 + 2/i 2 + ^) - 144/3 2 /i 2 /i 4 

24 (4 + 24/3/i/i + 3^ 4 (2/3 2 - 5)) ~ ' (3 ' 16) 

, _n [/3/i 3 (4 + /i 4 ) - 6h (4 + ^ (2/? 2 - 3)) - 24^/i 2 /x] 
M ~ 2 (4 + 24/3 V + 3// 4 (2/3 2 - 5)) ' 

where ' denotes differentiation with respect to r and where we have defined a rescaled dif- 
fusivity parameter 

/9 = -, (3.17) 
a 

which will play the role of the sole parameter in this system. This is a quantity of mass 
dimension one and we will denote as /3 the numerical coefficient of the Planck mass in this 
ratio: 

M a \ 3 ( M n 



^{wjywj- (3i8) 

where M H is the mass scale associated with x. 

Stability of Bounce Let us now show under which conditions the bounce is stable. As 
we will demonstrate, there is a range of values of f3 for which a stable bounce can take place. 
First, let us evaluate h' at the bounce point: 

_ ^ 2 (/ + 4) (15^-4) 
k ° ~ 96^ • (3 - 19) 

This is positive provided that // 4 > 4/15 and the no-ghost condition may also be satisfied 

A) = i + ^ 4 (V-2) >0 > ( 3 - 2 °) 

from which two conditions arise 

(3 2 > - and /U 4 > , (3.21) 
or f<\ and ^ < . (3.22) 

As we will show, the models with /3 2 > 5/2 do not bounce stably (see the condition arising 
from positivity of sound speed, Eq. (3.27)), therefore we have the requirement that for a 
ghost-free bounce 

k <lf< v^w (3 ' 23) 
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In this model, /i 4 = /x 4 , is the minimal value of /n 4 beneath which the energy for the scalar 
becomes positive at the bounce and would need to be compensated by a negative external 
energy density: this is the boundary of the dynamically inaccessible region, 

£o = f Q/4-l)=0 =► A = \. (3-24) 

The lowest value of /3 2 is such that the maximum ghost-free ^ 4 as given by the inequality 
Eq. (3.23) is also the minimum /^ 4 as given by Eq. (3.24), i.e. 

/3 2 > /? g 2 = I ■ (3.25) 

The upper bound on /3 2 comes from considering the sound speed at the bounce. For 
the hot G-bounce model, the positivity condition (3.11) reduces to 

fj, 8 + f - 3/3 4 ^) - 2/i 4 (9 - 4/3 2 ) + 4 > . (3.26) 

Again, the condition that there be at least one trajectory which has a stable bounce is 
equivalent to asking that the bounce be stable when the external energy density vanishes, 
i.e. // 4 = (j,q. Substituting this into Eq. (3.26) gives a condition on 0: 

2 <Pl = \ (j + Vu) . (3.27) 

For exceeding this upper bound, the sound speed at the bounce is always negative. In 
summary, stable bounces occur when 

1.29 ~ g < |0| <0 C ~ 2.07. (3.28) 

This condition is independent of the equation of state of the external fluid, provided that 
w > — 1. This is the case since we have derived this condition by looking at the limiting 
trajectory where the external energy density vanishes. 

One may then ask whether it is possible to bounce while keeping the perturbations 
subluminal, c 2 ^ 1. The answer is of course yes, since for close to C there are very few 
trajectories where sound speed squared is at all positive and for all of them it is close to 
zero. Let us be more precise: one needs to look at the full sound-speed expression, Eq. (3.5), 
which in the case of the hot G-bounce reduces to 

2 16 + 8^ 4 (40 2 -9)+ / u 8 (45 + 160 2 - 120 4 ) 



(4 + 3/i 4 (20 2 - 5)) 2 



(3.29) 



at the bounce point. This can be solved by requiring that a trajectory which bounces when 
the sound speed is exactly equal to the speed of light is stable to obtain a new limiting 
value, 14 



ft^ 49 ^ 1 , (3.30) 

Then trajectories which bounce stably with subluminal sound speed exist provided that 

1.64 ~ X < |0| <0 C ~ 2.07. (3.31) 

One should note that the sound speed typically increases following the bounce and can easily 
become superluminal. We should stress that this does not result in any causal paradoxes 
[37] even though it does signify that the UV completion of this model would not be Lorentz 
invariant [54]. 

14 The remaining solutions are not relevant since they fail the other stability tests. 
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Further Evolution Having established that it is possible to bounce stably in the hot G- 
bounce model, let us now turn to analysing the phase space numerically. We refer the reader 
to figure 1, where we have plotted the interesting part of the phase space. The striking 
feature of the hot G-bounce model is that there exist trajectories which after bouncing reach 
a maximum Hubble parameter and turn around to evolve toward the origin of the phase 
space — Minkowski spacetime. 

For any value of h there are limiting values of fj,, for which the external energy density 
vanishes (the boundaries of the red regions in figure 1) beyond which the system cannot 
evolve. The dynamically accessible region then is carved out by the inequality 

3h 2 > o£ = -/i 2 - + 3/i/?/x 2 . (3.32) 
2 8 

In the vicinity of the origin, /i, /i <C 1 so (3.32) reduces to 

h>-^= + 0(^), (3.33) 

where we have picked the trajectories in the lower-right quadrant of figure 1. 

The first type of trajectory we will consider is one corresponding to a universe with 
no accompanying external fluid: the evolution proceeds exactly along the boundary of the 
dynamically inaccessible pink region in this case. Taking the solution along the boundary 
close to the origin as h — — m/v6j we obtain as an approximation for the system (3.16) 

h' = -3h 2 + 0(h 4 ) , (3.34) 

We can solve this approximate system to obtain 

h — , r — > oo , (3.35) 

3r 

2 

which is an expanding universe comprising a fluid with a constant equation of state wx = 1- 
This is unsurprising, since at late times, when // < 1, only the leading-order term X in 
the Lagrangian function K is relevant. Thus the evolution is just that of a kinetic-energy- 
dominated canonical scalar field. 

In order to analyse the approach to the origin of the other trajectories, we turn to 
the first integral of the equations of motion, Eq. (3.4). In the vicinity of the origin the 
leading-order contribution to the first integral is 

I.t^. (3.36) 

This implies that close to the origin 

H ~ -(3Z) 3/4 /i 3/2 , (3.37) 
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Figure 1: Plot of the phase space for the hot G-bounce model (3.12) in the presence of 
external radiation fluid (w = 1/3), fi = ^/a<p, h = y/aH, j3 = x/a = 1.5. The solution is 
fully under control only in the blue region: pink regions are dynamically inaccessible, while 
the white regions suffer from instabilities (e.g. negative sound speed squared, < 0, or 
ghosts, D < 0.) In red, we have marked an example trajectory: it enters the healthy region 
from one where cj? < 0, bounces, reaches a maximum Hubble parameter and then proceeds 
toward the origin. The equation of state of the scalar approaches wx = 1, so eventually 
it redshifts away, leaving only the accompanying radiation fluid. In blue, we have marked 
a nearby trajectory: here the Hubble parameter does not turn around, but the trajectory 
proceeds to a pressure singularity (D = 0) where H diverges (but both the scale factor and 
Hubble parameter remain finite). Only a narrow range of bouncing trajectories arrive safely 
in the late-time radiation-domination era; the rest end in singularities. 



since this is the only solution which conserves the first integral as the origin is approached. 
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Again, approximating the system (3.16) on such a solution, we obtain 




ti = -2h 2 - ^-l 3 / 2 h 3 + O^) , (3.38) 

yl = y|(-^) 5/3 + O(^ 10/3 ) . (3.39) 

Eq. (3.39) can be solved easily to obtain 

/3X\ 3 / 4 

H ~ - ( — J t~ 3/2 , r -> oo . (3.40) 

We can also solve Eq. (3.38) provided we change the time variable to the scale factor, a. We 
then obtain the solution 

h=^Ca 2 -\J^-j , a^oo. (3.41) 

with C some constant of integration. By bringing this to the more user-friendly form of the 
Friedmann equation, we obtain 

1 V3X 3 

h '*tt + cw> a ^°°' (3 ' 42) 

which represents a radiation-domination phase with a small correction coming from a fluid 
with u; = 1. 

We can also calculate the properties for the perturbations of the scalar in the vicinity 
of the origin on the above background solution. We obtain 

D ~ 1 - 6/3(3Z) 3/4 /i 5/2 , h->0, (3.43) 

ci^i+m^f 1 ^' 2 . 

This is again what we would have naively expected assuming that all the higher-order terms 
in K and G become irrelevant and we are effectively dealing with a kinetic-energy-dominated 
canonical scalar field in the presence of radiation. 

Perturbations and Possible Inflationary Stage In this paper, we have not dealt with 
the perturbation spectrum and how it could match up to the observational limits. Models 
featuring an NEC-violating stage in the early universe presented to date produce the wrong 
spectrum unless furnished with an additional scalar field which acts as a curvaton [76-78]. 
We discuss these models further in § 4.1. 

In the case of the hot G-bounce model presented here, we also have some stumbling 
blocks. Firstly, the change in the scale factor during the collapsing pre-bounce stage is 
actually very small: only a few percent when the moment of crossing c 2 = is taken as the 
comparison point. This is the case for the hot G-bounce model and all the models presented 
in figures 2 and 3. 15 A stable collapsing universe per se is not difficult to construct — it 
has always been the bounce itself that was problematic and we have solved this problem 



15 The pressure singularities, D = 0, occur at a finite scale factor and Hubble parameter. The only divergent 
quantities are H and sound speed. 
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here. However, we have been unable to graft a long-lived collapsing phase onto our bouncing 
trajectories in a simple way. We remain hopeful that given enough effort this will turn out 
to be possible. 

Even without the collapsing phase one could imagine producing the spectrum of per- 
turbations in a way akin to tachyacoustic inflation [83-87]: As the trajectory approaches the 
pressure singularity, the sound speed becomes very high and then decreases as the system 
moves away. A decreasing sound speed provides for an acoustic horizon that shrinks (in 
comoving coordinates) and under the right circumstances is able to solve the horizon prob- 
lem even for a non-inflating background. A sufficiently large causal horizon could be set up 
during evolution along a trajectory which comes very close to the singularity. 16 However, 
in the hot G-bounce model, trajectories which pass close to the pressure singularity have 
large and rapidly- varying slow roll parameter H j H 2 and sound speed which would give a 
naive two-point function that is not scale-invariant. At this point, one would need to ensure 
that the standard Bunch-Davies vacuum is even established and to what extent the standard 
calculation of perturbation freeze-out carries through. This may not be the case for such a 
rapidly changing acoustic horizon [88]. Whether it would be possible to generate an appro- 
priate spectrum through such a mechanism warrants further investigation and we leave it to 
future work. 

For the moment, as a simple fix, we would like to point out that nothing prevents us 
from adding in a positive cosmological constant to the hot G-bounce model. This has the 
effect of making the final state of the evolution a de Sitter rather than a Minkowski space- 
time. In figure 2, we show that bouncing trajectories of the type that we have discussed 
above are in fact still present in the phase space, provided that the cosmological constant is 
not excessively large. 

We could then imagine that this cosmological constant is in fact a result of the existence 
of a very flat potential for the kinetically braided scalar 0, flat enough to be effectively 
constant until significantly after the bounce. We can then replicate a standard inflationary 
mechanism following the bounce. Just as in the standard inflationary scenario, a flat enough 
potential will give the correct perturbation spectrum following reheating. One may then 
ask whether it would be possible to have a sufficiently long stage of inflation to explain the 
perturbation spectrum at large scales, but retain some of the perturbations generated during 
the bounce at sufficiently large scales; for possible observational consequences see e.g. [89]. 

As such, the hot G-bounce model would not provide an alternative for inflation, but 
would be a possible extension of the inflationary mechanism into the past. We do not, 
of course, have full control of the trajectory at all times into the past: all the bouncing 
trajectories originate in regions with gradient instabilities. This is similar to the situation 
with inflation where the trajectories are also incomplete to the past, implying the existence 
of an initial singularity [6]. However, in this model, the energy scale of the crossover from 
regions of instability is in fact arbitrarily low (controlled by the mass scale in a), perhaps 
allowing us to avoid having to deal with quantum gravity. 

3.2 Other Shift-Symmetric G-bounce Models 

In figures 2 and 3 we have provided the phase-space diagrams for a number of shift-symmetric 
G-bounce models, demonstrating the wide range of trajectories possible in this class. In 
all the diagrams, only the blue regions are under full control: we have marked as pink 

16 One would have to worry about remaining within the weak-gravity regime and away from strong-coupling. 
Picking the right mass scale for a to achieve this may not be possible for trajectories too close to the singularity. 
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the dynamically inaccessible regions and as white the regions plagued by instabilities. The 
supplied bestiary of models is by no means exhaustive: the bounces are simple to achieve 
and the stable violation of the NEC opens up a whole set of alternative expansion histories 
for the early universe. 

Hierarchy of Sufficient Conditions As we have argued above, a stable bounce is realised 
when all the inequalities (3.6), (3.8), (3.11) and (2.11) hold. Now one can ask how simple 
is it to satisfy all these inequalities? Are our examples above, see figures (1), (2) and (3), 
very special / fine-tuned? Here we provide simple sufficient conditions for a healthy bounce. 
Namely the following conditions guarantee a healthy bounce: 

• kk^x > 0; 

• matter is normal w ^ — 1. Apart from this requirement, the models are indifferent as 
to the equation of state of external matter; 

• the following hierarchy is satisfied at the bounce 

113 

K > mK m > - (itik) 2 > > -m 2 K mm > - - (m K ) 2 . (3.44) 

Indeed the first inequality from the left guarantees (3.6) and K + w (K — mK jm ) > 0. 
Further the last two inequalities yield D > therefore we get (3.9), i.e. Hq > 0. Now let us 
see whether one can have a positive sound speed. Using mKK m = 2Xkkx > we obtain 

m/o (m) > -3m 2 k 4 + 2 (imK, 

then using the second and third inequalities from (3.44) we obtain 

m/o (m) > -3m 2 /t 4 — 3 ^2mK, m — (mz-t) 2 ^ k 2 + &K 2 mK. m = . 

This hierarchy (3.44) of inequalities is sufficient but not necessary for a stable bounce. How- 
ever, having two free functions k and K it is easy to chose them to satisfy (3.44) for some 
range of m. In that case this system will bounce in a healthy fashion in this chosen range 
of m for any type of external matter with w ^ — 1. In particular it seems to be possible to 
choose G and K to satisfy the hierarchy for all m. Note that it is also easy to construct a 
theory which allows for a stable bounce but violates the hierarchy, see e.g. the system on the 
top panel of figure (3) where the violation is manifest. These sufficient conditions are rather 
intuitive and can help for a future engineering of bouncing cosmologies. 
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Figure 2: A selection of phase-space portraits for shift-symmetric G-bounce models with 
external matter. The classical evolution presented here is only under control in the blue 
regions: the pink colouring represents dynamically inaccessible regions and the white regions 
have unstable perturbations (either negative sound speed squared or ghosts). 
Top Panel: Hot G-bounce model with the X 2 term restored: K = X+X 2 /2-X 3 , n = 2f3X, 
(3 = 1.5 in the presence of an accompanying radiation fluid. Restoring the X 2 term does 
not substantially change the phase space presented in figure 1, provided its coefficient is 
somewhat smaller than one. 

Bottom Panel: Hot G-bounce model with an additional positive cosmological constant: 
K = —A + X — X 3 , k = 2(3X, A = 0.05, j3 = 1.5, accompanying radiation fluid. Despite 
the addition of A, bouncing trajectories still exist, but now proceed to a late-time de Sitter 
attractor, providing an inflationary stage following the bounce. 
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Figure 3: A selection of phase-space portraits for shift-symmetric G-bounce models with 
external matter. Color coding as in figure 2. 

Top Panel: G-bounce with destabilised Minkowski: K = —X — X 3 , k = 2f3X, f3 = 1.7 in 
the presence of stiff matter, w = 1. In this type of models, the trajectories generically begin 
in a pressure singularity and then cross c 2 = (red) after bouncing, or vice- versa (blue). 
These models contain de Sitter attractors/repellers similar to those described in Ref. [73] 
and trajectories end or begin there. 

Bottom Panel: bounce & recollapse model: K = —A + X — X 3 , k = 2/3X + 7X 2 , negative 
A = —0.05, /3 = 1.7, 7 = —0.04 with radiation, w = 1/3. The red trajectory begins in a 
singularity at large H. The expansion slows down until the universe begins to recollapse. 
The crunching is prevented, however, and the universe bounces. The trajectory then crosses 
c 2 = eventually evolving to an unstable de-Sitter attractor. The blue trajectory is the time 
reverse: first a bounce, then a recollapse. 
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4 Bounces in Models with Negligible External Matter 



Models in which the external matter can be considered negligible in the vicinity of the 
bounce provide another category of bouncing cosmologies with a significantly simplified, 
two-dimensional phase space. We will discuss such models in this section and specialise to 
the discussion of the so-called conformal Galileon model [68, 76, 77]. 

The coordinates for the phase space in the class are provided by the pair (4>,ra). How- 
ever, the Friedmann equation (2.8) can be used, at least in principle locally, to eliminate the 
field value 4>, by substituting it with a function of m and H, 

(j) = ^{m,H). (4.1) 

This gives a phase space similar to the one obtained in section §3 with coordinates (ra,H), 
with the advantage again that the bounce position in the phase space is very explicit. On 
the other hand, the Lagrangian may have complicated dependence on (j) and therefore the 
solution Eq. (4.1) may have many branches, and it may not even be possible to obtain a 
closed form version of it. Nonetheless, in the example considered below, such a procedure 
does bear fruit. 

The evolution in this phase space is described by 

—2>n(2H — Km) — 2mn r k + 2K ( k 
m = ^ >* ± , (4.2) 

k (K j, — ran ^) — n {ran m + ran j, + 3kH ) 
H= ! ' 2D" ' • 

with D defined in Eq. (2.11) and any dependence on the field value (j) eliminated through 
Eq. (4.1). The above results are too general to proceed further. We restrict our attention to 
an example already known from the literature: the conformal Galileon. 

4.1 Conformal G-Bounce Model 

The action for the so-called conformal Galileon model is [68, 76] : 17 



/ 



S&= / d^v 73 ^ 



3 

(4.3) 



-2f 2 e 2<t> X + %rX 2 + %rXD^) 
A 3 A 3 



where A and / are constants with mass dimension one. We will consider this model as 
being minimally coupled to gravity. The conformal Galileon can be rewritten in terns of our 
notation as 

K = -f 2 e^m 2 + ^ra 4 , (4.4) 
K= 2 ^m 2 . (4.5) 

We proceed with the analysis of background cosmology in this model by applying the 
method described in §4: we can eliminate the scalar field value through the Friedmann 
equation, 

20 3/ 3 m 4 + 12ff/ 3 m 3 - 6A 3 H 2 
6 2/ 2 A 3 m 2 " 1 ' 

17 For further development see [77] and [78]. 
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This quantity must be strictly positive. This implies that wherever the above function is 
negative will be a region of the phase space which is dynamically inaccessible. 

Equations of motion for H and m can now be easily derived using their general equi- 
valents (4.2). It is helpful to rescale the time variable which enters both m and H, 18 



i) 



3/2 

r . 



We then obtain a rescaled chemical potential, fj,, and a rescaled Hubble parameter h, 

/ A \3/2 /A\ 3/2 
m = filj) and H = h(-\ , (4.7) 

and we will denote differentiation with respect to this new time variable with a prime, ( )' = 
d/dr. This procedure greatly simplifies the equations of motion and yields a set of equations 
which are free of parameters: 



, /i 3 [(I2h 2 + 5hfi + n 2 ) -fi 3 (6h + fi)] -2/i 2 (3/i + /i) 
" " M ^ (2^ + 1) + 2h 2 ' (48) 



ti 



^ [n 3 (l2h 2 + 16fyi + 3/i 2 ) - 8h 2 fi] - 12h A 
2 (// (2/i 2 + 1) + 2/i 2 ) 



This observation implies that the phase space topology is independent of the choice of the 
parameters / and A and only the physical scales corresponding to various trajectories will 
differ as /, A are changed. 

We are now able to write down the conditions for stability Eqs (2.11) and (2.12) in the 
conformal Galileon model in terms of the rescaled variables: 



D = ^ + + 6/x 4 > , (4.9) 

2 



3(H 2 D)~ 2 [I2h 4 - 32/iV + 4/iV (14/i 2 - l) + (4.10) 
+ I6h/j, 7 (1 - 2/x 2 ) + /i 8 (3 - 4/x 2 (2 + /?))] > . 

It can immediately be seen that D > at all points. Close to the origin, values of h are 
constrained by the requirement to keep Eq. (4.6) positive, which gives us h < ii 2 j\f2. This 
implies that for all the trajectories in the dynamically accessible region, D = at the origin. 
For all trajectories leaving or entering the origin d/i/d/x = and using l'Hopital's rule one 
can see that D — >• on those trajectories. 

Let us now deal with the scales present in this model. Firstly, the physical quantities 
(in Planck units) are related to the rescaled quantities defined in Eqs (4.7) through relations 
of the type 

H 2 =(^]\ 2 , H=(^]\'. (4.11) 



J J \f. 

This implies that introducing a hierarchy A <C / will reduce the physical scales in the problem, 

f 2 
z pi . 



allowing larger regions of phase space to retain curvatures significantly below order M| 



18 As opposed to the rescaled variables in the hot G-bounce model defined in Eqs (3.15), here both h and fi 
remain dimensioned (and are presented in Planck units). 
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Secondly, as discussed in Ref. [76, Eq. (17)], there exists a strong-coupling scale, which 
in Minkowski can be estimated to be 



Altrong ~ A 2 e 2 * . (4.12) 

This implies that in the part of phase space immediately surrounding the dynamically in- 
accessible region (the boundary of which is the surface e 2 ^ = 0, see Eq. (4.6)), the strong- 
coupling scale is always very low and approaches zero. This means that it is the trajectories 
which approach the boundary of the dynamically inaccessible region that will cross into a 
strong-coupling regime for the perturbations and not necessarily just a trajectory which has a 
large value of h. It should be noted, however, that a proper calculation of the strong-coupling 
scale taking into account the background solution has not been performed, thus one should 
not at this stage treat the strong-coupling problems calculated in this way as definite. 

Conformal Bounce Conditions Specialising the result Eq. (4.8) to the bounce point, 
h = 0, we obtain 

^ = J^> . (4.13) 

Since in the conformal Galileon model this quantity is always positive, the bounce point can 
only be crossed from collapse to expansion, if at all. Any trajectory that bounces will not 
proceed to recollapse, even if it may approach h = 0. 

The stability at the bounce point is determined by the sign of the quantity Eq. (4.10), 
which reduces to the condition 

4^ + 8/x 2 , - 3 < . (4.14) 

This can be solved, giving 

< |Mo| < ~y ~ 1 « 0.568 . (4.15) 
Substituting these limits into Eq. (4.13) gives us the constraint that during stable bounces 

0<h' < 17 « 0.095 , (4.16) 

or, equivalently in physical units 



A 



3 



< Hq < 10- 1 (jj Mpi . (4.17) 

The above shows that obtaining a stable bounce within the conformal Galileon model is in 
principle possible, as was noted in [76] and demonstrated in detail in [78]. 

We now turn to a global analysis of the phase space in order to understand the past 
and future of the trajectories which exhibit a bounce. We refer the reader to figure 4 and its 
caption for a description of the dynamics. 

Evolution of Trajectories Firstly, we comment on the behaviour of the evolution close 
to the origin of the phase plot which is a highly degenerate singular point. From physical 
considerations it is clear that the origin is a fixed point. There are no other finite fixed 
points for this system. In particular, in the dynamically allowed region close to the origin, 
\h\ < // 2 /\/2 + C(// 3 ). We can then find an equation describing trajectories in the vicinity 
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of the origin, h = h*((J,), by solving order by order for the coefficients of a Taylor series 
approximating the solution of 

dh _ /u 3 [fi 3 (I2h 2 + 16/i/i + 3^ 2 ) - 8h 2 fi] - 12/i 4 

d^ = 2fi (/i 3 [(12/i 2 + 5hfi + ii?) - n 3 (6/i + /x)] - 2h 2 (3h + ^)) ' (4 ' 18) 

This allows us to obtain two solutions. The first is the equation for the boundaries of the 
dynamically inaccessible region which, as has already been stated, are also solutions. Taking 
h< 0: 

/ l , b = -^ + ^-^ + ... = / , 3 -^7TTv. (4.19) 

On this trajectory, the system of equations (4.8) near the origin reduces to 

M ' = ^L + £?(//), (4.20) 

If we induce the Minkowski space to collapse, then the universe evolves on the approx- 
imate solution 

a ~ \= — - , r -> -co (4.21) 

h ~ — , 
2r ' 

which corresponds to a collapsing radiation-domination era. On this trajectory the perturb- 
ations have the following properties: 

D = 6/x 2 + 0(fi 3 ) , |U — >■ (4.22) 

c b 2 = ^ - £ + 0(m 3 ) • (4.23) 

A trajectory such as this one has been marked as blue on figure 4 and was studied in 
[78] where the above evolution corresponds to the initial pre-bounce stage of the solution. 

The second solution to (4.18) corresponds to a separatrix between the bouncing traject- 
ories in the left-upper quadrant of figure 4 and those trajectories which never collapse, but 
begin and end in a Big Rip and fill the top-right quadrant of the figure: 

h *z = Y + To^ + 3V + ° { ^ ] ■ (424) 

This trajectory is in fact the Galilean Genesis trajectory of [76] . We can again approximate 
the system of equations (4.8) on this trajectory 

M ' = /i 2 + 0(/z 4 ), (4.25) 
h> = \^ + O(^) . 

If we start slightly away from Minkowski along the Genesis trajectory, then we obtain an 
evolving solution which is always expanding. Eqs (4.25) can be solved: 

fi ~ — - , t — y — oo , (4.26) 

r 

h ~ 

2r3 • 
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Figure 4: Phase portrait for the conformal Galileon model, Eq. (4.3) in rescaled coordinates 
h = (f / A) 3 / 2 H , fi = (//A) 3 / 2 (/>. The main figure shows the stable region where /i > 0. The 
inset depicts the time-reversed region with fi < 0. The solution is under control fully only in 
the light blue regions: pink corresponds to dynamically inaccessible regions, white — to negat- 
ive sound speed squared. Yellow and orange are regions where curvature is transplanckian for 
(//A) 3 = 1, 2, respectively. The blue line is a typical healthy bouncing trajectory (presented 
in [78]): it originates from a region where the theory is strongly coupled, but the background 
solution evolves as a collapsing radiation-dominated cosmology; the universe then bounces 
in a healthy region and then the trajectory very rapidly crosses into the region where c 2 < 
and the classical solution should not be trusted. The red trajectory is the Galilean Gen- 
esis trajectory [76]: it begins in the vicinity of the Minkowski origin; the universe is always 
expanding and eventually the trajectory crosses the line c s = around h = 0.6; depend- 
ing on the choice of parameters this happens either before or after the curvatures become 
transplanckian. Both the trajectories merge to an attractor which evolves toward a Big Rip 
singularity. In the inset in green, we have marked a trajectory time-reversed with respect to 
the blue discussed above: this one begins in a (collapsing) Big Rip singularity, at some point 
crosses into a region of positive sound speed squared, bounces and then proceeds to expand 
in a radiation-domination- like phase which is also strongly coupled. 
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On these trajectories in the vicinity of the origin then 



D = 3/i 2 + 0(fi A ) 







(4.27) 




(4.28) 



This matches the asymptotics obtained in [76]. We have plotted this trajectory as the red 
line in figure 4. 

If we evolve the Galilean Genesis trajectory forward, it will cross into the region with 
negative sound speed squared. This occurs at h c ~ 0.6. Depending on the value of the 
Lagrangian parameters /, A, this may happen before or after gravity becomes strong, but is 
inevitable for all the trajectories. Indeed, all trajectories which have bounced in the past 
also eventually approach and cross the boundary c 2 = 0. 19 

On figure 4 it is hard to see that the transition through this boundary c 2 = occurs. 
This curve is given by the solution of the polynomial equation F (/i, h) = (/i 2 Z)) 2 c 2 = 0. If 
some trajectories (of nonzero measure) do not cross c 2 = then there should be an interval 
of h where these trajectories approach this boundary with dfi/dh = —Fh/F „, i.e. system of 
equations F (/x, h) = and dfi/dh = —F^/F^ should have a continuum of real solutions. As 
it is easy to check this is not the case and this polynomial system has 22 isolated complex 
roots. Thus all trajectories do evolve across c 2 = 0, because this curve is neither a trajectory 
nor a singularity of the system (4.8). 

Importantly, this Galilean Genesis trajectory does not evolve toward the region where 
the perturbations are strongly coupled as naively defined in Eq. (4.12), but in fact evolves 
away. Thus it is either strong gravity or gradient instabilities that will signify that the 
solution is failing, unless the effective field theory leaves its regime of validity. 20 In fact at 
the beginning of the Galilean Genesis trajectory, where h = and \x = 0, the naive strong 
coupling scale Agtrongj given by Eq. (4.12), is exactly zero. One would expect that this 
indicates that the scalar field is infinitely strongly coupled there. 

Moreover, when one integrates the Galilean Genesis trajectory numerically, one finds that 
the system effectively spends all of its time in the vicinity of the origin. Once h begins to 
pick up, the evolution is extremely fast and the Big Rip is reached in a short time. This 
means the time during which the system is in the region where h is significantly different 
from zero and yet there are no gradient instabilities is actually very brief. Hence the scale 
factor increases only by approximately 50% between the beginning of the trajectory in the 
vicinity of the origin and the trajectory's leaving the stable region, i.e. only about 0.4 e- 
folds of expansion are under control for the perturbations. In the Galilean Genesis scenario 
effectively all the expansion occurs in the region where gradient instabilities are present in 
the degree of freedom driving the expansion. 

Also in the bouncing trajectories of the type (4.19) the amount of expansion between 
the bounce and the crossing into c 2 < is quite small. However, here all the perturbations 
are set up during the collapsing phase, which can effectively be arbitrarily long. 

One should be concerned, however, since the naive strong coupling scale (4.12) is in 
fact extremely low in the whole region where the contraction mimics radiation-domination 
and the perturbations are generated. The proper strong-coupling scale for a cosmological 

19 In fact all the bouncing trajectories enter the region where perturbations are out of control at an even 
smaller value of h, h < 0.6 

20 See the discussion in [72] regarding the rather subtle question of the appropriate cut-off for the effective 
field theory. 
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background should be calculated in order to judge this with any certainly. This is outside of 
the scope of this work. Despite the fact that the observable perturbations are driven by a 
curvaton, a separate scalar field, the predictivity of the calculation for curvaton perturbations 
is also uncertain. The curvaton by design is coupled explicitly to the conformal galileon, thus 
if the galileon perturbations are strongly coupled, the curvaton perturbations will also not 
be under control. 

5 Conclusions 

In this paper, we have demonstrated that stable bouncing cosmologies are generic and simple 
to achieve in models containing non-canonical scalar fields with Kinetic Gravity Braiding. We 
have constructed models that can realise a transition from contraction to expansion in a flat 
Friedmann universe stably, without leaving the weak-gravity regime. The ingredients can — 
but do not have to — include external matter, which will generically blue shift as the universe 
contracts. 

The conditions for the bounce to be healthy are that the perturbations remain ghost 
free and that the evolution is free of gradient instabilities. Indeed, we have derived a set of 
sufficient conditions on the form of the Lagrangian of the theory which will ensure that the 
evolution at the turnaround point itself is stable. As we have shown, constructing Lagrangians 
which contain healthy bouncing trajectories is not difficult. 

On the other hand, we have not succeeded in constructing a model where the whole 
expansion history, including the remote past and remote future, is under control. We have 
found that, generically, all trajectories which bounce usually have some kind of pathology. 
A typical bouncing trajectory begins or ends in a Big Rip or pressure singularity. At best, 
the trajectories cross the line of vanishing sound speed, corresponding to a singularity in the 
acoustic metric, which can appear in non-canonical theories [37] . Since quantum fluctuations 
are normalised by the sound speed, they will grow without bound as the trajectory approaches 
the singularity. However, this is a singularity in the scalar sector, the energy scale of which 
is essentially arbitrary. Thus, these trajectories avoid being transplackian and as such do not 
necessarily involve strong-gravitational effects. 

We have succeeded in finding a category of models which, while starting from such an 
acoustic singularity, bounce stably and then evolve to a healthy future, the precise nature of 
which depends on the accompanying fluid. The late-time solutions can correspond, for ex- 
ample, to a hot Big-Bang — when the external fluid is radiation — or to a standard inflationary 
stage — if the scalar field is furnished with a very flat potential which breaks shift-symmetry 
weakly. 

Despite this success, the bouncing models we have proposed do not resolve the problem 
of the initial singularity of the universe, similarly to inflation [6]. One could say that we 
have shifted back the beginning of the universe's clock to a collapsing, pre-inflationary stage. 
However, we remain optimistic that one should be able to eventually resolve this issue. The 
models we have presented give an idea of what a solution could look like: a model with a 
fixed point inside a stable region could have trajectories which are healthy throughout their 
complete history and never evolve to a strong-gravity regime. The ease with which one can 
construct models featuring stable bounces gives hope that some additional complication (for 
example, interactions of the braided scalar and other degrees of freedom) will result in a 
realistic bouncing model of the universe. 
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